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Abstract 

The coefficients of the nonlinear terms in a modified Altarelli-Parisi evolution 
equation with parton recombination are determined in the leading logarithmic (Q^) 
approximation. The results are valid in the whole x region and contain the transla- 
tion GG — > qq, which is inhibited in the double leading logarithmic approximation. 
The comparisons of the new evolution equation with the Gribov-Levin-Ryskin equa- 
tion are presented. 
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1 Introduction 



The inclusion of parton recombination to the proton structure function through a mod- 
ified AltareUi-Parisi equation is an interesting subject involving the test of perturbative 
QCD and the study of new effects in small x physics. A traditional tool in this research is 
the so-called Gribov-Levin-Ryskin (GLR) equation, presented in two pioneering papers: 
one is the idea of shadowing arising from gluon recombination, which was proposed by 
Gribov, Levin, Ryskin [1] based on the AGK (Abramovsky, Gribov, Kancheli) cutting 
rules [2] in the double leading logarithmic approximation (DLLA); the other is a pertur- 
bative calculation of the recombination probabilities in the DLLA by Mueller and Qiu [3] , 
which enables the GLR equation to be applied phenomenologically. 

Unfortunately, a series of questions in the GLR equation destroyed the gluon recombi- 
nation effects. In our previous work [4], one of us (WZ) has proposed that the application 
of the AGK cutting rules is unreasonable in the GLR equation since it sums up the dia- 
grams in which the cut lines break the important correlation among the initial gluons. For 
this reason, a new evolution equation including parton recombination was established in 
the leading logarithmic {Q^) approximation (LL((5^)A) using time ordered perturbation 
theory (TOPT) instead of the AGK cutting rules [4] . This equation for parton distribu- 
tions in the proton can be generally written as 




dG{xB,Q^) 



,GG-^GG 



(xi,X2,X3,X4, A) 



5{xi -\- X2 — x^ — X/^)[5{x^ — xb) + 5{x4^ — X B)\dxidx2dx^dx^d/^ 
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-2{^f I ^2, XI + A, X2 - A, Q^) E P?''^''''{X,, X3, X4, A) 

-ftC^ Jx\>xb j 

5(a;i + a;2 — a;3 — a;4)[5(a;3 — xb) + 5(a;4 — a;s)](ia;i(ia;2(ia;3(ia;4(iA, (la) 

and 

= (^)' / G(2)(^i, ^2, a;i + A, 0:2 - A, Q2) E i^^^^^^(a;i, 0:2, 0:3, x,, A) 

6{xi + a;2 — a;3 — a;4)[5(a;3 — xb) + S{x4 — x B)]dxidx2dx3dx4dA 
-2(^)2 / G^'\x^,X2, X, + A, a;2 - A, Q^) ^ /^^^--^^(xi, 0:2, 0:3, 0:4, A) 

-flC^ Jxi>Xb j 

5(xi + ^2 — X3 — X4) [5(3:3 ~ 3:5) + 5{x4 — XB)]dxidx2dx3dx4dA, (16) 

where we temporarily neglect the linear terms in the Altarelli-Parisi equation and G^^^ 
is the correlation function for four-gluon lines. The new equation provides the following 
physical picture for the gluon recombination in a QCD evolution process: the two-parton- 
to-two-parton (2 2) amplitudes (Figs, la-lb) lead to positive (antiscreening) effects 
and the interference amplitudes between the one-parton-to-two-parton (1 — > 2) and the 
three-parton-to two-parton (3 — > 2) (Figs. Ic-lf) amphtudes lead to negative (screening) 
effects, respectively. The TOPT-analysis shows that the above two kinds of amplitudes 
correspond to the same recombination function J2i Pi but in different kinematic ranges. 
On the other hand, the complete corrections of the gluon recombination to the evolution 
equation should include the contributions from the virtual diagrams corresponding to 
S{xi — Xb) and S{x2 — xb), however, they cancel each other in the proton [4]. 

The quantitative predictions of Eq. (1) depend on the recombination functions. 
Mueller and Qiu [3] have derived similar transition probabilities J2i PMQ,i in the DLLA, 
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where only the gluon ladder graphs are kept at small x. We should remember that the 
AGK cutting rules are apphcable only in diagrams which consist of gluon ladders [1]. 
Therefore, the DLLA is a necessary approximation for using the AGK cutting rules in 
the GLR equation. On the other hand, in the Mueller-Qiu approach [3], the cut vertex 
technique [5] was used to calculate the transition probabilities, which are entangled in a 
complex cut diagram with four-parton propagators [3]. Therefore, the DLLA is also a 
convenient approximation in Ref. [3]. 

However, we shall point out that the DLLA leads to a difficulty in the GLR equation: 
the transition of gluon— > quarks is suppressed in the DLLA-manner. Although the spe- 
cial techniques are used to include the corrections of gluon recombination to the quark 
distributions in Ref. [3], however, it seems inconsistent in the theory as we will show in 
Sect. 5. 

The purpose of this paper is to derive the recombination functions of Eq. (1) in 
the LL((5^) approximation. We can generalize the recombination function to the whole 
X region and even include the processes with quarks, since Eq. (1) is not restricted 
by the AGK cutting rules [4]. Furthermore, we used TOPT instead of the cut vertex 
method in the derivation of Eq. (1) [4]. Thus, we can separately calculate the simple 
two-parton-to-two-parton (2 — > 2) processes in the recombination function. Following the 
above mentioned derivation, we shall complete a new modified Altarelli-Parisi evolution 
equation with parton recombination in the LL((5^) approximation, which is valid in the 
whole X region. Momentum conservation is naturally restored due to the coexistence of 
the shadowing- and antishadowing-effects in the new evolution equation. We also find 
that the DLLA is a bad approximation even for the recombination function of gluons in 
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the small-a; region. 

The outhne of the paper is as follows. In Sect. 2 we shall derive a set of complete 
recombination functions. An undetermined factor in Eq. (1) is the correlation function 
G^'^\ which relates to the nonperturbative structure of the proton. We use a simple 
assumption to model G^"^^ in Sect. 3. Combining the results of the above two sections, we 
establish a modified Altarelli-Parisi equation including parton recombination in Sect. 4. 
In Sect. 5 we compare our new equation with the GLR equation. 

2 Recombination functions 

The recombination function in Eq. (1) is factorized due to the application of TOPT 
in [4]. All processes of two-parton-to-two-parton (2 — > 2) type can be precisely calcu- 
lated using the standard perturbative QCD, except that the energies of partons are not 
conserved at the vertex connecting with the probe. 

At first, we calculate the recombination function for the two-gluon-to-two-gluon pro- 
cess. The momenta of all initial and final partons are on-shell and they are parametrized 
as (see Fig. 2), 



Pi = [xip, Q, xip] ; p2 = [x2P, Q, X2p] , 



p'l = [x'lP, Q, x[p\ ; P2 = W2P, Q, X2P] , 



X3P+ 7r^^^±^^3P 

2xsp 



X4P+ , -l±,X4P 

IX4P 



(2) 



We take the physical axial gauge and the light-hke vector n fixes the gauge as n- A = 0, 
A being the gluon field. The corresponding recombination function is defined by [4] 
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:[M(piP2 - ki')[M{p[p', ki')ru 



,1.0 d l 



Ej,, + Ej,, — Eh — Ell 



^2Ek' {2T:f2Ev 

The index i in Eq. (3) implies the t-, u-, s-channels and their interference terms. We 
begin with the i-channel graph (Fig. 3). The contribution of this process to the invariant 

amplitude is 



\MtM; 



X3=Xb 



„4 fi2 -1 -1 KiKj illiim 



where {i,j),{l,m),{p,q) and (r, s) are the space-indices corresponding to (a,/?), 
(Cv) (0)X)) respectively. We need to distinguish the probing place, for example, the 
transfer momenta and Ir are determined by two down- vertices in Fig. 3 where the 
probing place is X3 = x^, i.e.. 



k 
Ir 



{X4, - X2)p + 



I' 



2X4P 



{X4 - X2)p 



{xi - x'^)p ■ 



2xip 



, -/x, (X4 - X'^)p 



(5) 



The algebra in Eq. (4) can be straightforwardly derived without any approximation 
in computer. The contributions from the w-channel and the interference channels can be 
similarly obtained by using the interchanges of the corresponding momenta. 



Now we turn to discuss the s-channel as shown in Fig. 4. We should note that the 
massless partons with the parallel momenta can go on-mass-shell simultaneously in the 
coUinear case and the coUinear singularity may arise in the s-channel since = (pi+P2)^ = 
(p'l +^2)^ — 0. Fortunately, we now have an useful tool [6,7] to pick up the short-distance 
contributions in the propagator with coUinear divergence: we use the following special 
propagators (for quarks) 



and (for gluons) 



to replace the normal Feynman propagators, respectively. Using the definitions 



n^ = (nO,nx,#) = -^(1,0^,1); n'^ = -^(1, 0^, -1), (8) 

or equivalently 



n'^ = (n+,n-,7Tx) = (1,0,0J; n'^ = (0, 1, 0^), (9) 

we have = 0, = 0, and n ■ n = 1. The contribution of the s-channel (Fig. 4) to the 
invariant amplitude can then be safely calculated as. 



\MsM*\ 



9^ 



CI 



(^) 

4 m -ri-n^ 



uiuj 111 Jim 
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where we have used the following replacement in the calculations: 



(11) 



Obviously, the contributions given by the s-channel and corresponding interference 
channels integrate to zero if the correlation function G^'^\xi,X2,x[,X2,Q'^) is symmetric 
under xi X2 and x[ <-> X2, since these contributions contain the factor {x2 — xi) or 
{x'2 — x[). Thus, we can forget the s-channel and its interference terms in GG — > GG. 
The contributions of Feynman diagram with four-gluon vertex are neglected because they 
are less singular. The total contributions of GG —>■ GG are given in Appendix A. 

The recombination function Y,i P^'^'^'^'^ has poles at xi = and X2 — Q (see Eq. 
(A.l)). However, Eq. (l.a) is really infrared safe. One can simply check this conclusion 
as follows. For example, Xi = implies that A = since only twist-4 or twist-2 amplitudes 
contribute to an unpolarized amplitude. Summing over the contributions of the right-hand 
side of Eq. (l.a) at the poles, one can find that 




>GG — >GG 



{xi = 0, X2, A = 0, X3, X4, A) 



S{Xi + X2 — X2 — X4)[d{x3 — Xb) + 5{X4^ — XB)]dX2dXzdX4^ 



+ 



1>x\>xb 



G^^\x,, a;2 = 0, A = 0, Q^) ^ i^^^^^^(a;i, 0:2 = 0, A = 0, 0:3, x^, A) 



5{xi + X2 — x^ — X4)[5{xz — Xb) + 5{x^ — XB)]dxidxzdxi 



2 / G^^H 



J \>X\>Xb 



xi,X2 = 0,A = 0,Q2)^p.< 



,GG^GG 



(xi, a;2 = 0, A = 0, xs, x^, A) 



6{xi + X2 — — X4)[6{x3 — xb) + 5(2:4 — XB)]dxidx3dx^ = 0, (12) 

where the two positive terms are symmetric under xi ^ X2 and Xi 7^ in the negative 
interference term. Therefore, Eq. (1) is infrared safe. 

Furthermore, unhke the GLR equation, the momentum is conserved in Eq. (1), i.e.. 



djQdXBXBG{xB,Q^) _ ^ , . 

and 



dJo dxBXBq{xB,Q^) ^ , 

dlnQ2 ^ ^ 

Now let us return to consider the recombination functions for GG qq. For example, 

the contribution of the i-channel process in Fig. 5 to the invariant amplitude is 



X3=Xb 



4 iV(iV"= — Ij X2X2 l± 

The final results for X^f*^^^^ a-re hsted in Appendix A, where the contributions of the 
s-channel and corresponding interference channels vanish due to the same symmetry as 
in Eii^^^^^^. 

We need only consider the gluon recombination at small x, since gluons dominate the 
small-x behavior of the parton distributions in proton. However, in principle, our equation 
can include the recombination of quark-quark, quark-antiquark and quark-gluon in the 
LLA. For this end, we list all recombination functions in Appendix B. We will discuss the 
properties of our results in Sect. 5. 
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3 Parton correlation functions 

Equation (1) includes the gluon correlation function G^'^\xi,X2; x'l, x'2, Q^). In general, 
the parton density is a concept defined at twist-2. The parton correlation function is a 
generalization of the parton density beyond the leading twist and it has not yet been 
determined both in theory and experiment. For comparison with the GLR equation, we 
shall use a toy model as in Ref . [3] . First of all, we assume that 

G^'^\xi, X2, xi + A,X2- A, g^) 

^ G^^\xi,X2, xi + A,X2- A, Q^)S{xi - X2)d{A). (16) 

One can understand Eq. (16) as follows. Parton recombination is expected to be signifi- 
cant only when the actual spatial overlap between the fusing partons becomes sufficiently 
large in the interaction time Tint of probe-target both in the transverse- and longitudinal- 
directions. This means that (i) the transverse area occupied by partons should become 
compared to the effective transverse area; (ii) only two partons with the same value of 
rapidity y = In^ and same impact parameter can overlap their wave functions in the 
longitudinal direction during Tint, since the initial partons are parallel moving along the 
2;-direction. Therefore, the fusion between the gluons with different values of y is sup- 
pressed. 

We consider that the different branches of the parton cascade evolve independently; 
this is valid in the large limit [8]. Thus, G^'^\x,Q^) ~ G\x,Q^), where G{x,Q^) 
is the usual gluon density. However, the density G{x, Q^) can not be normalized since 
G{x, Q^) ~ 00 when x ^ 0. A more reasonable approach is to use the gluon density in 
rapidity space xG{x, Q^) — dna/dy instead of G{x, Q"^). The normalization condition is 
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due to momentum conservation, where the sum is over all parton flavors. Therefore, we 
can assume the two-gluon distribution per unit area in proton to be 



{xi + X2)G^^\xi,X2, A, = ^^[x,G{x,, Q')][x2G{x2, g')]5(xi - X2)5{A), (18) 

where the first factor is from a normalization factor of the parton correlation function [3,6] 
and the second factor is determined by a simple geometric model as in Ref . [3] . Although 
Eq. (18) is a phenomenological model, however, its x-dependence can be checked by 
experiment. 

4 Modified Altarelli-Parisi equation 

In this section we derive a modified Altarelli-Parisi equation with gluon recombination. 
Inputting Eq. (18) with the recombination functions to Eq. (1) and adding the linear 
terms corresponding to the Altarelli-Parisi equation, we have 

dxBG{xB,Q'^) 



dlnQ2 

r ^^P^-<^ix„xn)x,G{^^,Q') 

TT JxB Xi Xi 



^l^i^'f j^'J dXiXBX,G\xuQ'')^P°°^''ixuXB), {19.a) 

where the recombination functions Eqs. (A.1)-(A.2) become simple after integral over 
X2,A., xs and X4 under the assumption (18): 
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i 

_ Sal C\ (2x1 - xb) [-\S<qxbx\ - 64a;ia;| + \S2x\x\ + mx\ + 16x|) 

which is vahd in the whole x region. The second term on the right-hand side of Eq. 
(19. a) is positive due to the contributions of Figs, la-lb and is called as antishadowing 
effect. The negative term in Eq. (19. a) arises from the shadowing corrections from the 
interference diagrams Figs, l.c-l.f. 

We have a similar result for GG — * gg: 

dxBS{xB-,Q'^) 
= ^ r ^-^^P^^\x,,xb)xM^uQ') 

TT JxB Xi Xi 

r^'dx,XBX^G\x,,Q')Y:P['''~'\x,,XB), (19.C) 
IbTT^ RQ JxB i 

where the recombination function with the assumption (18) is 



j:Pr-\xr,XB) 

i 

Tf (2xi-XBf {4^x1 + 5x1- 6xiXb) 
N{N^ - 1) 

^ NTf {2xi - XBf{4:xl - QxiXb + 4:xj] 



[19.d) 



(N^ - 1)2 xl 
We will discuss the properties of our new modified Altarelli-Parisi equation (19) in the 
next section. 
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5 Discussions 



In order to compare Eq. (19) with the GLR equation, we begin with a review of 
the DLLA. As definited by the AltareUi-Parisi equation, the DLLA means that in each 
order in ctg, one keeps only the ln((5^///^) ln(l/x) factor in the solutions of the evolution 
equation, or equivalently, only the terms having 1/z = xi/xb factor in the splitting 
function are necessary to generate larger logarithms in (Fig. 6). We know that the 
splitting functions PAp'^i^) -^1p*^(^) have the lowest 2;-power (~ l/z), while PXp'^i^) 
and Pap'^{z) vanish in the DLLA. Therefore, the DLLA diagram consists of the gluon 
ladders, since any transition of gluon— > quark breaks the ladder-structure. Similarly, 
we also find (see the Appendices A and B) that only P^*^^^ , pm^G pqG^G g-^^ 
the leading contributions in the DLLA since they have the factor 1/z — {xi + X2)/xb- 
Thus, we can conclude that any transition of G — > g or GG — > g is suppressed in the 
DLLA-approach. 

Now let us remember the GLR equation, which has following form in [3]: 



dxBG{xB,Q^) 
dlnQ2 






and 



dXBS{xB,Q'^) 

dlnQ^ 
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-0.17 



160i?2g2 



asm 



j 



xb Xi Xi 



dXi Xb j-,GG^qq 



'xiH{xi,Q^), (20.6) 



with 



dlng2 




(20.c) 



and 



MQ 



(z) = -2z + 15z^ - 30z^ + 18/. 



(20.d) 



As we have emphasized, the DLLA inhibits any transition of G — > g or GG q. 
Therefore, Ref. [3] takes some special treatments to reahze the above mentioned transition. 
For example, the last term on the right-hand side of Eq. (20.b) with Eq. (20. c) comes from 
a contribution of the next leading logarithms in and is inconsistent with the DLLA. On 
the other hand, 5{xi —xb) is inserted by hand in the derivation of the second term of Eq. 
(20. b). However, 5{xi—xb) implies a cut line through an initial gluon line with momentum 
Xi and is really the contribution of the virtual diagram to {d/d\nQ^)xBG{xB,Q^)', this 
contribution should therefore be canceled in the proton [4] . Even if we forget the above 
mentioned two questions in Eq. (20), we still cannot prevent the linear term in Eq. (20) 
from mixing the quark line with the gluon ladder at each Q?. In consequence, either the 
evolution from to + AQ^ will stop in the DLLA-manner due to the inhibition of 
gluon— >quark, or the DLLA as well as the AGK cutting rules, which regard the gluon 
ladder as the pomeron [1] , can not be applied in the GLR equation. 

Unlike Eq. (20), our new equation (19) includes the transition of gluon— >quark in the 
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whole X region. Another important difference between Eqs. (19) and (20) is that the 
positive antishadowing effects are separated from the negative shadowing effects in Eq. 
(19) because they have different kinematical domains in x. 

We emphasize that the DLLA is an unreasonable approximation even in Eq. (19. a). 
In fact, we take the DLLA and keep the leading term ~ 1/z in Eq. (19. a). In this case, 
this equation can be simplified as 

dXBG{xB,Q'^) 

dlnQ2 

So;., 



TT JxB 



Jxb/2 xi HQ JxB xi 

where the upper-limit of the integral in the negative interfering terms is 1/2 due to the 
restriction X1+X2 < 1. There is such gluon density, for example, xG{x, Q'^) ~ x^''{l — xy, 
where the net recombination effects disappear due to the balance of the shadowing- and 
antishadowing-effects. The calculation shows that Xc — 1 and 0.5 for Eq. (19. a) and 
Eq. (21), respectively. Therefore, we can find that the DLLA obviously distorts the 
recombination effects even in the gluon evolution equation at small x. 

Finally, we estimate the size of the gluon recombination effect. On can use 

nonlinear terms 

W = — > a,, 22 

linear terms 

to determine the kinematic range, where the nonlinear effects can not be neglected in Eq. 
(19). The approximation solutions of Eq. (22) are 
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xG{x, Q^) > 



{RQf 



far Eq. (19.a), 



2.4 



{RQr 



far Eq. (19.6). 



(23) 



- 1.13 



The HERA data show that xG(x, Q^) ~ 2 - 10 for = (1 - 10) GeV^ and x = 10"=^. 
Thus, we can expect that the gluon recombination effects wiU be appeared in the proton 
structure function in the range < 10 GeV'^ and x< 10~^ if << 1 fm (i.e., the gluon 
distribution in the proton has the "hot spots" -structure). 

In summary, the coefficients of the nonhnear terms in a modified AltareUi-Parisi evo- 
lution equation with gluon recombination are determined in the leading logarithmic (Q^) 
approximation. The results are valid in the whole x region. The comparisons of the new 
evolution equation with the GLR equation are presented. We expect that the gluon re- 
combination effects of the new equation can be appeared in the proton structure function 
in the HERA-domain of small x and low Q"^ if the "hot spots" structure exists in the 
proton. 
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Appendix A 

We give the recombination functions in Eq. (1) for the case of Xi — x'l and X2 — X2 
follows. 

i 

■' i 

= 2«^ ^2 ^ 1 (^1 + a;2 - a;s)(255a;?a;^ + 255a;ta;^ + 146a;ia;^ + UQx\x2 + 270a;?a;^ 
+Uxlx% - l?,2xlxB - l'i&xlx% + 186a;^a;| - 132x^xs - 136x^a;| + 186a;^a;| + 64x^x| 
+58xi + 58^2 - 2Ad>xlx2X% + SlOx^XaXs + ^blxixlx\ - 2ASxixlx\ + ^^lxlx2x\ 

—AQ1x\x\xb — 2^?)x\x2Xb — AQ1x\x\xb — 288X1X2X5 + 64xiX2X^)/[8xbX^X2(xi + X2)^]. 

{A.l) 

One can get the general result for xi 7^ x'^ and X2 7^ x'2 if we replace xi and X2 by 
^Jxix'i and \jx2x'2 in Eq. (A.l). 

^ J^'^'^"^^(xi,X2,Xb) 
i 

= / 51 Pi^^~""^{^l: ^2, 2^3, X4)5(xi + X2 - X3 - X4)(5(X3 - XB)dXzdXi 

i 

+ / X]-^j'^*^^^^(^i'^2,a;3,a^4)(^(a;i +X2 - X3 - X4)5(x4 - XB)dxzdxi 

Xi + X2 — xbY{Ax\ + 2XiX2 + X2 - 2x2Xb + 2x\ — AxiXb) 



^2al\2 



N{m - 1] 



{xi + ^2)^X2 



^ (xi + X2 - xb)^(4x| + 2xiX2 + xf- 2xiXb + 2x1 ~ 4x2Xb) 

(Xi + X2)^Xi 
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NTf {xi+X2-XBf{^xl-^XiXB-^X2XB+xl + 2x2X-l + xl) \ 

(iV2_i)2 {xr + X2Yx2Xr j ' ^ - > 

where we don't distinguish quark and antiquark. In the above mentioned equations, we 
have not includes the contributions of the s-channel and corresponding interference terms, 
since they integrate to zero. 
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Appendix B 



For simplicity, we present the recombination functions of quark-quark and quark-gluon 
at Xi — X2 and x[ = x'^^ where the contributions of the s-channel are included due to the 
asymmetry of the two-parton distribution under xi <-> X2. 



. y x^ 



04: X-^ 



Y,P"{xi xb) = 1 (2^1 ~ XB){Qxj - xiXb + 12xbxI - xl) 



V- ^ N _ 7 {2xi - xb){10xi - xb) 



(S.4) 



V- „qG->G/ X 1 (2a;i - XBf{3Qxl + 62a;? - 39xbXi) 

l^Pi [xi,xb) = — 4 . {B.5) 

^ 144 Xbx\ 



^ pgg^G/ N _ -1 (22^1 - xb){^x\ - 25x1x1 - ixlxB + 18x1x1 - 9x%) 
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Figure Captions 

Fig. 1 The diagrams contributing to a new modified Altarelli-Parisi equation with gluon 
recombination. Here the shaded part implies the correlation of gluons at short-distance 
and "x" means the probing place. 

Fig. 2 The two-gluon-to-two-gluon subprocess in recombination function. 
Fig. 3 A t-channel diagram for GG — > G. 
Fig. 4 A s-channel diagram for GG — > G. 
Fig. 5 A t-channel diagram for GG — > q. 

Fig. 6 A dominant splitting process in the DLLA, where ^ Xi leads to a strong order 
in X. 
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